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O Abstract 

^ ,Anion-radical form of the oxygen centers O is predicted at the DFT level for small silver oxide particles having the AgO sto- 
ichiometry. Model clusters (AgO) n appear to be ferromagnetic with appreciable spin density at the oxygen centers. In contrast 
^**^ jto these clusters, the Ag20 model cluster have no unpaired electrons in the ground state. The increased O/Ag ratio in the oxide 
■particles is proved to be responsible for the spin density at oxygen centers. 
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1. Introduction 



£h In the last decades silver oxide has been actively studied in 
-i— > numerous works due to its importance in various applications. 

^ 'In particular, silver oxide is considered to be used in batteries 

q i (Hi. molecular sensors J2] and bactericidal materials |H,|3]. The 
. i-H 'interaction of oxygen with metal silver surface leads to the for- 

K^.mation of silver oxide overlayers, which play a crucial role in 
^ (— i the silver-catalyzed ethylene epoxidation and the partial oxida- 

£^tion of methanol to formaldehyde Si. 

I— ~~ '■ Photoelectron spectrum for in-situ prepared AgO bulk sam- 
I ,ples contains a double-peak Ols signal with a remarkable 
^> chemical shift of 2.9 eV between the two components [7]. The 
CN| authors assigned the splitting of Ols peak via the co-existence 
of two nonequivalent oxygen in - 1 and -2 oxidation states in sil- 
ver oxide structure. Later, it was also remarked that silver oxide 
thin films with high oxygen content contain the same oxygen 
forms ii. This is in agreement with earlier DFT calculations 
'for the AgO oxide in which allowed the authors to suggest ap- 
pearance of "holes" in electronic structure of O 2 oxygen |@]. 
■ The silver-catalyzed ethylene epoxidation reaction is com- 
. £h . mon ly considered to be determined by the competition between 
two oxygen forms with different Ols electron binding ener- 
5^ gies called electrophilic and nucleophilic oxygen species lioll . 
The electrophilic oxygen is suggested to be responsible for 
the epoxidation route of ethylene oxidation, while the nucle- 
ophilic oxygen governs the combustion 

[11]. 

For a long time 

the nature of electrophilic oxygen has been debated in terms of 
atomic B12I1 and molecular [13] forms. Although, in a number 
of works the nucleophilic and electrophilic oxygen species are 
proposed to be two atomically adsorbed forms with different 
charge states, both forms were usually considered as the O 2- 
species with closed shells. The first model of radical oxygen, 



o 



to the best of our knowledge, belongs to Carter and Goddard 
lfl4l [I5I1 . who introduced so-called surface atomic oxyradical 
oxygen (SAO) O in addition to the standard oxide oxygen O 2 
on base of a cluster Ag30. Later, the increase of electrophilic 
atomic oxygen concentration was experimentally shown to ac- 
celerate the ethylene epoxidation reaction, while increasing the 
concentration of nucleophilic oxygen slows it down B16I1 . 

The origin of atomic electrophilic oxygen on silver catalyst 
surface and its electronic state are still not completely under- 
stood. The electrophilic oxygen is known to be formed on the 
silver surface at saturation with oxygen accompanied by the for- 
mation of silver oxide overlayers [5] and has Ols binding en- 
ergy about 2 eV higher than that for the nucleophilic oxygen. 
This oxygen also is characterized by a lower effective charge 
about 1 a.u. 
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In our previous paper II18I1 closed and open shell states of 
oxygen adsorbed on silver were simulated based on the Ag20 
molecule to model nucleophilic and electrophilic species, re- 
spectively. Unlike Carter-Goddard SAO in Ag30 model hav- 
ing single unpaired electron at the oxygen center, the open- 
shell Ag20 was predicted to be of the singlet biradical type 
having two electrons with antiparallel spins separated between 
oxygen and silver. The increase of the Ols energy and disap- 
pearance of the pre-edge feature of Ox-edge XAS spectrum for 
the electrophilic oxygen were in a perfect agreement with the 
experimental observations. Since the open-shell oxygen form 
emerged in the excited state of the system, the question arises 
for which oxo silver systems the atomic electrophilic oxygen 
appears in the ground state. 

It is well-known that nanoscale transition metal oxides pos- 
sesses unusual structural, electronic, magnetic and catalytic 
properties, which are very different from the same properties 
of bulk samples lfl9ll . Nanostructured silver oxide is not an 
exception revealing extreme photoactivity [ 20] . However, re- 
liable experimental data on the atomic structure of silver oxide 



Preprint submitted to Chemical Physics Letters 



December 5, 2012 



clusters are currently lacking. As a rule, quantum chemical cal- 
culations up to now have been focused on the study of silver 
oxide clusters with regular stoichiometry Ag2<3 12111 or oxygen 
adsorption on silver clusters 11221 l23l 12411 .Silver oxide clusters 
with the stoichiometry AgO has not yet been studied theoreti- 
cally to the best of our knowledge. 

Nanostructured silver oxide AgO can have magnetic proper- 
ties due to the appearence of anion-radical oxygen in the struc- 
ture which is an obvious result of a deficiency of the ls(Ag) 
electrons to form the 2 ~ centers. In contrast to the interac- 
tion between silver and oxygen in stoichiometry Ag20, in AgO 
each oxygen atom is able to withdraw only one silver 5s elec- 
tron giving rise to the oxygen state 0~. Formation of the O" 
state in AgO is schematically illustrated in figure Q] Although 
mixed oxidation state of silver atoms might be also realized, it 
is not the case for small AgO clusters as will be shown below. 




Figure 1 : The scheme, which illustrates the appearance of oxygen species 2 ~ 
and CT due to the interaction between silver and oxygen with stoichiometry 
Ag20 and AgO, respectively. 



In case of silver oxide AgO unpaired electrons are to be local- 
ized on oxygen atoms, contrary to the majority of the magnetic 
transition metal oxide systems. Although magnetism of small 
AgO clusters has not been experimentally proven yet, isoelec- 
tronic copper oxide nanoparticles are known to exhibit ferro- 
magnetic or paramagnetic properties which is explained by the 
increased O to Cu stoichiometric ratio II251 12611 . 



In the present work, the existence of anion-radical form of 
oxygen centers in the ground state is demonstrated for the 
model silver oxide clusters (AgO) n (n = 1-4, 6). The properties 
of oxygen in the clusters (AgO) n are compared with the prop- 
erties of oxygen in simple silver oxide clusters having various 
stoichiometries to show how an increase in the relative content 
of oxygen leads to an increase in spin density at oxo center. The 
electrophilic oxygen centers on silver responsible for the epoxi- 
dation is suggested to have anion-radical structure found for the 
small (AgO) n clusters. 



2. Methods 

All quantum chemical calculations in the present work were 
carried out within the density functional theory using the Am- 
sterdam Density Functional package i27l 281. The combi- 
nation of the exchange functional Perdew and Wang 1986 
(PW86x) I29II and Perdew and Wang 1991 correlation func- 
tional (PW91c) l30ll has been applied using polarized triple- 
zeta basis set of Slater-type orbitals (TZP). To take into account 
relativistic effects ZORA 13111 was additionally used for silver 
atoms. 

Geometry optimizations of all clusters were performed, start- 
ing from structures having reasonable Ag-O bond lengths. Op- 
timizations were carried out without restriction of symmetry. 
The optimized geometries appear to have positive vibrational 
frequencies which means that these structures correspond to lo- 
cal minimums on the potential energy surface. For all struc- 
tures the states with different spin projections S ; were checked 
in spin-polarized calculations. 

3. Results and discussions 

The smallest silver oxide cluster AgO has calculated bond 
length of 2.06 A, which is in a good agreement with the experi- 
mental value of 2 A [32]. The structure of (AgO)2 is a rhombus 
with the Ag-Ag distance of 2.91 A, the O-O distance of 3.27 
A, and the Ag-0 distance 2.19 A. Experimental bond lengths 
for Ag2 and O2 molecules are 2.69 A and 1.21 A, respectively. 
Thus, both the Ag-Ag and O-O bonds in (AgO)2 are weaker 
than their homonuclear counterparts. The (AgO)3 cluster has 
planar ring-like structure. In thus cluster the Ag-Ag bond is 
significantly smaller than the O-O bond, and, consequently, the 
structure of this cluster is not a perfect hexagon. The average 
Ag-O bond length in (AgO)3 is 2.14 A, being close enough to 
that value for (AgO)2- The distorted cubic structure is real- 
ized for the cluster (AgO)4 containing the cluster (AgO)2 as a 
building block. The average Ag-O bond length is 2.31 A, be- 
ing markedly larger than that in smaller clusters. The Ag-Ag 
and O-O distances of 3.31 A and 3.22 A, respectively, become 
close to each other. The structure of (AgO)6 is simply an ex- 
tension of the (AgO)4 geometry with the addition of a (AgO)2 
unit with the formation of distorted prism. Although the Ag-O 
bond lengths vary between 2.16 and 2.34 A, the average Ag-Ag 
and O-O distances remain close to the corresponding ones in 
(AgO)zt. Optimized geometries of the clusters (AgO) n (n = 1-4, 
6) in the ground states obtained in the calculations are shown in 
figure |2] 

Most metal oxides are predominantly ionic in its structural 
features B19I1 . So, calculated small silver oxide clusters (AgO) n 
have the same structural motifs as some other well-known ionic 



clusters, for instance, NaCl, MgO, CsF, etc ||33i|34|]. Ag-O bond 
lengths in (AgO) n with n=l-3 are in magnitude the same as in 
different bulk silver oxides 1I35I1 . Ag-O distances in (AgO)4 and 
(AgO)e are similar to the Ag-O bond lengths in silver oxide thin 
films or oxygen adsorbed on silver surface 13611 . 

The stability of clusters (AgO)„ with n > 6 was estimated by 
the value of the atomization energy (E a ) defined as 
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Figure 2: The structure of the clusters (AgO)„ (n = 1-4, 6) as revealed at the PW86-PW91/ZORA-TZP level. Arrows point out the localization of unpaired electrons 
in the ground state (ferromagntic) of clusters. 



E a (n)=E(Ag)+E(0)--E n 
n 

where E(Ag), E(O) are the total energies of Ag and O atoms, 
and E n is the total energy of the cluster (AgO) n . 

In this definition, E a is positive for clusters that are stable 
with respect to the complete atomization. The value of E a in- 
creases monotonically with increasing size n, implying the pos- 
sibility of larger clusters formation (table [TJ. For instance, in 
our work [37] dependence of atomization energy of (ZnO) n 
(n=2-9) of the size also increased monotonically and larger 
clusters have been observed experimentally l38ll . 

Considered (AgO) n clusters appear to have an electron struc- 
ture in which unpaired electrons are localized almost exclu- 
sively on oxygen centers. Spins on oxygen are coupled through 
non-magnetic silver centers to form a ferromagnetic structure. 
Such spin ordering differs from that for magnetic oxides in 
which spins are localized on metal centers while the oxygen 
centers serve as non-magnetic bridges as shown in particu- 



lar in previous computational studies on the magnetic prop- 
erties of transition metal oxides clusters |34, 39[ 4(J ■ Figure 
[3] shows spin density distribution for the (AgO) n clusters in 
ground states. 

The low-spin states of the clusters (AgO) n under consider- 
ation are less stable than the high spin states. These states 
might be assigned to antiferromagnetic ordering of spins on 
oxygen centers. Such magnetic structures are given by broken- 
symmetry solutions revealed by the distinct large value of the 
spin contamination, i.e. the excess of the mean value for the 
spin-squared operator over that for the given value of the spin 
projection S^Sj+l) 114 lL 14211 . The spin contamination in DFT 
corresponds to the negative spin density or the spin polariza- 
tion as was shown by Wang, Becke and Smyth 64311 . In DFT, 
the status of spin contamination is quite different from that of 
the Hartree-Fock theory (where it is an indicator of the spin- 
symmetry purity) since the Kohn-Sham (KS) determinant is not 
a correct wave function and so could not be contaminated in the 
sense of Hartree-Fock theory. The spin contamination calcu- 
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Figure 3: Spin density distribution for the clusters (AgO) n (n = 1 - 4, 6) in ground states as revealed at the PW86-PW91/ZORA-TZP level 



lated for the KS determinant indicates in what extent the KS or- 
bitals are polarized. It is worthwhile noting that KS orbitals are 
always polarized (unrestricted) for open-shell molecular sys- 
tems as the effective potentials for the one-electron equations 
differ between a and /3 electrons as was pointed out by Pople 
and co-authors 844(1 . The spin contamination in DFT might be 
explored as the effective number of spatially separated electrons 
with antiparallel spins: if the spin contamination is about 1.0 
then there is one such pair, 2.0 - two pairs, etc 114 111 . The ex- 
ception is the Agf,06 cluster in the quintet state, in which there 
are four unpaired electrons having parallel spins almost equally 
delocalized over six oxygen atoms. (AgO)2 cluster has not a 
singlet state with a rhombic atomic geometry. 

In the case of the cluster (AgO)6 ferromagnetic and antifer- 
romagnetic states are almost degenerate. In the different spin 
states system has virtually the same geometry with a slightly 
different bond lengths. This indicates that the unpaired elec- 
trons in this system interact weakly with each other. One may 
assume that with an increase in cluster size, the low-spin anti- 
ferromagnetic state is to become energetically more favourable 
than high-spin ferromagnetic state. 

Thus, in the clusters with stoichiometry AgO in both the 
high-spin and low-spin states the oxygen centers appear to be 
of the anion-radical type O . With the increase in the cluster 



size the charge and spin density on oxygen is reduced reflecting 
gradual shift of the electron density to silver. 

One can expect that bonding mechanism in the (AgO) n clus- 
ters is predominantly ionic and is determined by electron-pair 
formed between 5 s electron of Ag and 2p electron of O. How- 
ever, the anion-radical oxygen has unpaired electron on itself 
implying that 3-electron bonding can be realized. In order 
to better understand bonding mechanism in the (AgO) n model 
systems bond energy analysis available in the ADF program 
l45l H3l is applied. The interaction energy AE; nt are decom- 
posed into three physically meaningful terms: 

A£ int = AV e i st + AEpauH + A£oi 

which correspond to the classical electrostatic interaction en- 
ergy, the Pauli repulsion energy and orbital interaction energy, 
respectively. The 3-electron bonding energy, which can be de- 
termined as AEp au ii + AE ;, contributes about 20-30% to the 
total bonding energy of (AgO) n systems, while the classical 
electrostatic interaction energy makes up 70-80% of the total 
bonding energy. It is the 3-electron bonding that is responsi- 
ble for the stabilization of the anion-radical oxygen in (AgO) n 
clusters. 

To compare the properties of oxygen in (AgO) n clusters 
with those of the clusters with other stoichiometries, the Ag20, 
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Table 1: The size of cluster n, the atomization energies E a , spin-projection Sj, mean value for the spin-squared operator (S 2 >, Mulliken charges p and spin densities 
p s for the ground (high-spin) states of clusters, as well as those for excited (low-spin) clusters. 



n 


1 


2 


3 


3 


4 


4 


4 


6 


6 


6 


6 




0.50 


1 


1.50 


0.5 


2 


1 





3 


2 


1 





<s 2 > 


0.76 


2.02 


3.77 


1.77 


6.03 


3.01 


1.99 


12.03 


6.13 


2.76 


1.82 


E a , eV 


4.227 


5.317 


5.814 


5.780 


5.855 


5.771 


5.802 


6.158 


6.146 


6.145 


6.144 


i* 


Mulliken charges p on atoms 


1 


0.547 


0.629 


0.689 


0.690 


0.656 


0.662 


0.653 


0.775 


0.869 


0.923 


0.888 


2 


-0.547 


0.629 


0.625 


0.618 


0.656 


0.691 


0.653 


0.648 


0.646 


0.633 


0.627 


3 




-0.629 


0.625 


0.624 


0.656 


0.663 


0.669 


0.649 


0.649 


0.636 


0.629 


4 




-0.629 


-0.662 


-0.680 


0.656 


0.668 


0.670 


0.648 


0.648 


0.609 


0.619 


5 






-0.661 


-0.671 


-0.656 


-0.683 


-0.645 


0.648 


0.648 


0.607 


0.620 


6 






-0.615 


-0.582 


-0.656 


-0.669 


-0.645 


0.786 


0.888 


0.928 


0.942 


7 










-0.656 


-0.662 


-0.678 


-0.700 


-0.731 


-0.743 


-0.746 


8 










-0.656 


-0.669 


-0.677 


-0.699 


-0.732 


-0.742 


-0.741 


9 
















-0.679 


-0.714 


-0.680 


-0.674 


10 
















-0.679 


-0.716 


-0.686 


-0.689 


11 
















-0.697 


-0.727 


-0.743 


-0.739 


12 
















-0.698 


-0.729 


-0.742 


-0.736 


i* 


Mulliken spin densities p s on atoms 


1 


-0.030 


0.065 


0.149 


0.151 


0.109 


-0.032 


-0.001 


0.299 


0.145 


0.005 


0.024 


2 


1.030 


0.065 


0.131 


0.057 


0.109 


0.117 


-0.002 


0.108 


0.117 


0.094 


0.041 


3 




0.934 


0.124 


0.051 


0.109 


0.018 


0.121 


0.109 


0.118 


0.097 


0.014 


4 




0.934 


0.878 


0.827 


0.109 


0.112 


-0.122 


0.108 


0.118 


0.063 


0.001 


5 






0.870 


0.830 


0.890 


0.239 


0.787 


0.108 


0.117 


0.063 


-0.014 


6 






0.848 


-0.917 


0.890 


0.865 


0.791 


0.311 


0.121 


0.006 


0.013 


7 










0.890 


-0.181 


0.870 


0.010 


0.582 


0.354 


0.066 


8 










0.890 


0.862 


-0.871 


-0.005 


0.585 


0.376 


-0.216 


9 
















0.893 


0.486 


-0.452 


0.779 


10 
















0.893 


0.477 


0.653 


-0.478 


11 
















0.797 


0.569 


0.361 


-0.072 


12 
















0.795 


0.563 


0.378 


-0.158 



*number of atom in cluster (see fig. [2} 



Ag3C>2 and AgC>2 systems were calculated. Their structure is 
shown in figure [4] Table [2] shows the values of the spin pro- 
jection S ; , mean value of spin-squared operator (S 2 ), Mulliken 
charges p and spin densities p s on atoms and energies of Ols 
Kohn-Sham orbitals for a and f3 sets. 




^Ag «0 



Figure 4: The structure of the clusters Ag20, Ag3C>2 AgC>2 as revealed at the 
PW86-PW91/ZORA-TZP level 

With the oxygen stoiciometry increasing in the sequence 
among Ag20, Ag3C>2, AgO and AgC>2 clusters the negative 
charge on oxygen decreases causing the drop of the Ols Kohn- 
Sham levels. This finding agrees well with the experimental 
data by Kaichev with co-authors who suggested that the charge 
on oxygen is responsible for the relative value of the binding 



energy of oxygen Ols electrons on the silver surface lfl7ll : the 
increase in the charge on oxygen decreases the binding energy 
of its Is electron. 

Oxygen in Ag20 cluster having zero spin density and a sig- 
nificant negative charge is similar in its properties to the oxy- 
gen in bulk Ag20 silver oxide with formal oxidation state of -2. 
The same properties are inherent to atomic nucleophilic oxy- 
gen species which are routinely registered during adsorption of 
oxygen on regular silver surface. 

In turn, anion-radical oxygen appearing in small (AgO) n 
clusters has significant spin density and smaller magnitude of 
negative charge. The binding energy of the Ols level of the 
anion-radical oxygen is significantly stabilized as compared to 
that of the oxygen in Ag20 cluster. All these oxygen parameters 
seem to correspond to the oxidation state of - 1 . 

It is O oxygen that was suggested to reveal itself in the pho- 
toelectron spectroscopy of bulk silver oxide AgO. At a high 
coverage oxygen may intiate a reconstruction of silver sur- 
face into AgO-like thin oxide film. It may cause the forma- 
tion of anion-radical oxygen center on surface which may gov- 
ern silver-catalyzed ethylene epoxidation. The latter sugges- 
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Table 2: Mulliken charge p and spin density p s on atoms and energies of Oi s Kohn-Sham orbitals (for a and p sets - in spin-polarized cases) of the clusters Ag20, 
Ag 3 2 ,AgOand Ag0 2 





S 2 


S, 


P 




E KS (0 ls ) (a;j3), eV 


Ag 


O 


Ag 


O 


Ag 2 








0.41 


-0.82 








-509.54 


Ag 3 2 


0.79 


0.5 


0.49 


-0.72 


-0.17 


0.52 


-509.85; -509.50 


AgO 


0.76 


1 


0.54 


-0.54 


-0.03 


1.03 


-511.42; -510.89 


AgO z 


3.76 


1.5 


0.77 


-0.38 


0.23 


1.38 


-513.52; -512.75 



tion is supported by the fact that the electron acceptors such as 
chlorine facilitate the ehylene epoxidation reaction (allowing 
the formation of electrophilic oxygen at low oxygen concentra- 
tions), while the electron donors such as cesium slow down the 
reaction [47]. 



4. Conclusions 

In this Letter, we performed density functional study on 
small silver oxide clusters (AgO) n (n=l-4, 6) to demonstrate 
the possibility of the formation of anion-radical oxygen partic- 
ipating in the formation of silver oxides. 

The results clearly indicate that in all considered clusters 
with stoichiometry of AgO the oxygen centers are radical-like. 
The spins on oxygen are ferromagnetically ordered for small 
clusters, while for the larger clusters antiferromagnetic state is 
likely to be more favourable. 

The spin density and charge on radical-like oxygen centers 
in the (AgO) n particles is larger than that in Ag 2 cluster. 

The anion-radical oxygen is suggested to have the same 
nature with atomic electrophilic oxygen species which govern 
silver-catalized ethylene epoxidation and the oxygen form in 
the oxide AgO. 
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